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The Hubbard model in bipartite lattice has been rigorously proved to have 77-paired state as its 
eigenstates in Refi. In this paper, this result is first generalized to triangular lattice, and then even 
arbitrary lattice structure. Then we also consider its generalization to higher-spin Hubbard models. 
The conditions under which these models have r;-paired states as its eigenstates were obtained. 



I. INTRODUCTION 

The fermionic model can possess off-diagonal long 
range order(ODLRO) and thus becomes superconduct- 
ing via particle pairing^'. The usual BCS state does have 
ODLRO via the mechanism of Cooper pairing, but it is 
not an eigenstate of a Hamiltonian with a local interact- 
ing term. Thus, it is of long-standing interests to find out 
models whose eigenstates possess ODLRO. One of the ex- 
act results about this has been obtained in Ref.^ , in which 
a local two-particle pairing operator called ?7-pairing was 
constructed and proved to be an eigen-operator of spin- 
1/2 Hubbard model in square lattice. Thus the states 
created by the 77 operator are exact eigenstates of the 
Hubbard model and display ODLRO. 

As one of the few rigorous results about the Hub- 
bard model in higher dimension, the 77 pairing has at- 
tracted a lot of theoretical interests since it was pro- 
posed. Although in the strong repulsive interaction case 
these eigenstates created by the 77 operator, which are all 
doubly-occupied states, have a large excitation energy, 
and therefore is not closely related to current discussion 
of high-Tc superconductivity of Cuprates in the content 
oit — J model, the rj operator has been found very help- 
ful and widely used in the study of attractive Hubbard 
model^^. Although in the majority of cases the 77-pairing 
state is an excited state, it has been found to be ground 
state for some extended Hubbard models, such as in pres- 
ence of nearest-neighbor interaction, bond-charge inter- 
action and so onSi. Besides, the properties of ?7-pairing 
operator have also been used to construct an SU{2) sym- 
metry group of the Hubbard model^. Together with 
another SU(2) group generated by particle-hole opera- 
tors, an 5'0(4) symmetry of the Hubbard model has been 
found, which unifies the spin order and superconductor 
pairing order^. 

Recently, rapid theoretical and experimental develop- 
ment has been made in the field of ultracold atoms, 
among which the successful application of optical lattice 
technique is a remarkable one^ The fermionic atoms 
confined in optical lattice can be described by (extended) 
Hubbard models. Using the technique of quantum optics, 
various lattice structures can be realizedii. It was also 
proposed that the tunnelling elements between nearest 
lattice sites can also be experimentally controlled so that 
the atoms in optical lattice can be effectively viewed as 
charged particles moving in a magnetic fieldiSi^. In this 



paper, we will try to generalize the ry-paring to the tri- 
angular lattice, and even to arbitrary lattice structures. 

On the other hand, the lattice atomic gases also pro- 
vide a testing grounds to study high-spin extension of 
Hubbard models, and the interaction between atoms is 
spin-dependent and can be tuned in a wide range due 
to the Feshbach resonance technique. Recently, Wu, Hu 
and Zhang discussed interacting spin-3/2 fermion atoms 
in optical lattice. They found that this model not only 
possesses a generic 5*0(5) symmetry formed by particle- 
hole operators, but also has an SU{2) symmetry con- 
structed by T] paring operators in the condition of some 
special interaction parametersii^ We will also discuss the 
77-pairing for more higher spin models in the following. 

We organize this report as follows. In the second sec- 
tion we will briefly review the keys to obtain the com- 
mutation relation [r'i'',H] — Crj'' for spin-1/2 Hubbard 
model in square lattice, where H is the Hamiltonian and 
C is a constant. In the third section we will general- 
ize it to Hubbard model in triangular lattice. Because 
the triangular lattice is not a bipartite lattice, the 77 op- 
erator is generally not an eigen-operator of the kinetic 
energy term. However, we find for a class of special 
hopping elements, a modified 77 operator can be eigen- 
operator. In the fourth section, we will discuss the in- 
teracting fermions with higher spins. We will answer the 
question that how to choose the interacting parameters 
so that an extended 77 operator can be an eigen-operator. 



II. THE 7] PAIRING 

Before discussing the generalization, we should first 
remind ourselves of some basic ideas of the 77 pair- 
ing discussed in Refi for the Hubbard model on two- 
dimensional square lattice, whose Hamiltonian is written 
as 

H = T+V = -ty^^{a]^ai^j+a],-aij+h.c.) + Uy^^n-(inii. 

ij i 

(1) 

A two-particle pairing operator of spin-1/2 fermions can 
be generally written as 

77 = > (7(fc)a -a 

k 



(2) 



it is an eigen-operator of the interaction energy term 
when (Ti'iyJTi') — (7i|y|77) is a constant for any \n) and 



|n'), for which (ri|?7'''|n') is nonvanishing. Here the states 
\n) are all eigenstates of the interaction term and they 
form a complete bases of the Hilbert space. Noticing that 
the interaction terms are all local, the pairing should also 
be a local one, which means that g{k) should be taken as 
a constant independent of k. 

In the momentum space the kinetic energy term can 
be written as 



f = 5^/(^)4 



Or 
k,ij k,a- 



kcr 



with 



f{k) — —2t{coskx + cos ky). 



(3) 



(4) 



In Refi it is noticed that the rj operator commutes with 
the kinetic energy term when q is taken as if because 



cos k + cos(7r — fc) = 0. 



(5) 



Thus in the coordinate space the rj pairing is a superpo- 
sition of local particle pairs 



V 






(6) 



where one essential point is the sign of the pair in one site 
should be different from those in its nearest-neighbors. 
Hence the rj operator thus defined satisfies the commu- 
tation relation 



r]^H - Ht]^ = -2U7]K 



(7) 



III. 



GENERALIZATION TO TRIANGULAR 
LATTICE 



We notice that the above discussion about rj pairing re- 
lies on whether the lattice structure is bipartite, because 
the lattice structure directly determines the explicit form 
of /(fc), which plays a crucial role in finding rj operator. 
Hence, for the Hubbard model in triangular lattice, the i] 
paring operator introduced in previous section is not gen- 
erally an eigen-operator. Then the question arouses that 
under which condition the model has a pairing operator 
as its eigen-operator. 

The Hamiltonian of the Hubbard model in a triangular 
lattice under consideration is written as 



H = j2*^jala,, + UY, 



ni-iUii 



(8) 



Here 9i are the phases defined in links, whereas only 

n 



closed loop 



are gauge invariant quantities, which rep- 
resent the flux number through the closed loops. In the 
momentum space, the kinetic energy term can be written 
as 

fk^ a/3 ^ 

cos {kx + Oi) + cos -^ + —ky + 02 




-I- cos 



We notice that when 91—62 — 9^ — 7r/2, the expression 
of f{k) turns out to be 

/(fc) = -4icosf^ + ^-^ 
■'^ ' V 2 2 4 

- cos ——k,, 



In this case, because 



/(fcx, ky) + f[2Tr- kx, -^{93 -92)-ky]=0 



(11) 



(12) 



the ry operator defined as 77^ = ^ a -a ^ - with q taken 

k 

as {2tt, -7^(6*3 — 92)) is an eigen-operator of this model. 

Therefore we have succeeded in obtaining a sufficient 
condition for the rj pairing states being eigenstates, which 
is 9i — 92 — 93 = 7r/2. A typical configuration is plotted 
in Fig^ This condition in fact means that there is a 7r/2 
fiux through each plaquette. 

Furthermore, we notice that the rj operator is a lo- 
cal operator, and the kinetic energy term only involves 
nearest-neighbors hopping, the assumption of transla- 
tional invariance is indeed not essential. Let iy = |t|e'^*^ , 
the condition can be further extended to 



E' 



-(modTr), 



(13) 



First focusing on the case with translational invariance, 
the kinetic energy term T can be written in a general 
form 

ia 



where C stands for any closed loop with odd numbers of 
links. This criterion can be easily verified in the coordi- 
nate space representation, and it is valid for any other 
lattice structure. One can see that the bipartite lattice 
is just a special case of this criterion, because any closed 
loop in a bipartite lattice must contain even numbers of 
links. 
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FIG. 1: Two typical configurations of flopping efements tij, 
wfiicfi satisfy tfie requirement Eg. 11311 . Tfie number 1 or i de- 
notes tfie flopping efements tij form one site to anotfier site. 
For tfie upper one, the Hamiftonian fias transfationaf invari- 
ance, and for the lower one it is not transfationaf invariant. 



IV. GENERALIZATION TO HIGH SPINS 
FERMIONS 



K,5J,,]+ = 0. 



(17) 



The first term describes the interaction between Tn/2 
state and — m/2 state, the second describes the interac- 
tion between ±7n/2 states and ±r7z'/2 states with dif- 
ferent m and to', and the last term describes the spin 
flipping processes in the channel of total Sz equalling to 
zero. The interacting between two spin-3/2 atoms dis- 
cussed in Refii^ is a special case of this form. One can 
find that the presence of the spin-flipping term is a con- 
sequence of spin-dependent interaction. 

It is illuminating to consider the commutation relation 
between alaj and these terms. For the first term, the 
commutation relation is 



[alal, ala|a|aj_ 



For the second term we have 



~a\a\. 



-ajajfejjfoo 



(18) 
(19) 



and for the spin flipping term 

[ajaj, 6J6|a|ax]_ = -h\h\a\a^ - h\h\a\ai + h\h\. (20) 

Hence, to satisfy the requirement rj'^V — Vrj^ ex ry^, the 
terms such as a!aJ6j.6cr should be cancelled with each 
other, which means that the rj operator should be taken 



In this section we will consider an interacting spin- 
^"""'^ fermions model, with 2n different spin states in 
each site. In this model, the two-particle pairing operator 
is generally not eigen-operator of the interaction term 
because of the enlargement of the Hilbert space. It is our 
interest to find out an extended f] operator and when the 
pairing operator is an eigen-operator of the interaction 
term V . 

We consider the extended Hubbard model as 



H^^tija\^aj„ + V, 



(14) 



ija 



where iy satisfy the criterion Ea. (|13|) . It can be shown 
that a better choice of the on-site interaction should be 
like 

V = V2{a\a^a\ai + b\b^b\bi + c\c]c\ci H ) 

+yQ^^{a]jaab\,b„> + ala„c\,c„i + b\b„c^^,Ca' H ) 



(7,(T' 



J J} 



J J 



+Vi{a\a\b^bi + h.c. + ajajc|C| + h.c. H ). (15) 

Here we denote a| and ai as fermionic annihilation op- 
erators for ±1/2 states, 6| and 6| for ±3/2 states and so 
on, i.e. 



[aa,a'^,\ + 



Sarr' 



(16) 



= V e*'"(at,at, + bi,bi, + ct,ct, +•••), (21) 



and the parameters in Eq. H15|l should satisfy the require- 
ment 



Vi = -2^0. 



(22) 



This relation recovers the result obtained in Ref.^'*. The 
commutation relation is therefore 



[v\V]^~{V2 + 2Vo)7i\ 



(23) 



The generalized 77-operator is also a superposition of local 
spin-singlet two-particle pairing. The state created by 
the operator has ODLRO for the same reason discussed 
in Ref-iand Refi^ 



V. CONCLUSION AND DISCUSSION 

In summary, we have discussed two generalizations of 
the 77-paring to two kinds of extended Hubbard mod- 
els. The original conclusions about the ry-pairing was 
restricted for bipartite lattice structure, in this paper we 
showed that for the Hubbard model defined on any lat- 
tice structure, a modified ?7-operator can be the eigen- 
operator if the hopping elements satisfying the criterion 



Ea. (|13|l . We also discussed the r^-pairing for higher-spin 
Hubbard models. We find that in order that the model 
has modified 77-paired states as its eigenstates, the in- 
teraction part should contain spin-flipping term and the 
parameters should satisfy the condition Ea. (|22() . 

As the original 77-paired stated, the generalized 77- 
paired states are also usually excited states. One thing of 
much concern is when these generalized 77-paired states 
become ground states. In the original paperi, it is shown 
the ry-paired state is metastable for attractive interac- 
tion. And later, various extended models which have the 
77-paired state as their ground states were proposed^^. 
Following these ideas, how to obtain models having the 



generalized ?7-paired states as the ground states is under 
further consideration. We hope that these efforts could 
shed light on the study of fermionic atoms in optical lat- 
tice and some solid-state materials such as NaCoO which 
have triangular lattice structurei^. 
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